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ABSTRACT 

Using the new capabilities of the Spitzer Space Telescope and extensive multiwavelength data from 
the Spitzer Infrared Nearby Galaxies Survey (SINGS), it is now possible to study the infrared prop- 
erties of star formation in nearby galaxies down to scales equivalent to large H ii regions. We are 
therefore able to determine what fraction of large, infrared- selected star-forming regions in normal 
galaxies are highly obscured and address how much of the star formation we miss by relying solely on 
the optical portion of the spectrum. Employing a new empirical method for deriving attenuations of 
infrared-selected star-forming regions we investigate the statistics of obscured star formation on 500 pc 
scales in a sample of 38 nearby galaxies. We find that the median attenuation is 1.4 magnitudes in 
Ha and that there is no evidence for a substantial sub-population of uniformly highly-obscured star- 
forming regions. The regions in the highly-obscured tail of the attenuation distribution {Ana ^ 3) 
make up only ~4% of the sample of nearly 1800 regions, though very embedded infrared sources 
on the much smaller scales and lower luminosities of compact and ultracompact H ii regions are al- 
most certainly present in greater numbers. The highly-obscured cases in our sample are generally 
the bright, central regions of galaxies with high overall attenuation but are not otherwise remarkable. 
We also find that a majority of the galaxies show decreasing radial trends in Ha attenuation. The 
small fraction of highly-obscured regions seen in this sample of normal, star- forming galaxies suggests 
that on 500 pc scales the timescale for significant dispersal or break up of nearby, optically-thick dust 
clouds is short relative to the lifetime of a typical star-forming region. 
Subject headings: galaxies: photometry - galaxies: ISM - stars: formation - H ii regions 
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1. INTRODUCTION 

Stars form enshrouded in their parent molecular 
clouds, but in time clusters of stars emerge from the ob- 
scuring dust and gas as H ii regions and disperse part 
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of this natal material. In the simplest case, the emis- 
sion that we observe from such systems should undergo 
a transition from being strongly attenuated in the opti- 
cal and ultraviolet (UV) to being less attenuated as the 
stars begin to break free of their birth clouds. 

In reality, the amount of obscuration observed in re- 
gions of star formation can vary widely depending on 
the environment. In nuclear starburst galaxies, V-band 
attenuations (Av)^"^ as high as 10-100 magnitudes are 
commonplace. The most well-studied of all, M82, has 
a central starburst with Ay ~25-30 (Rieke et al. 1980; 
Puxley 1991; Mattila & Meikle 2001), and the central 
starbursts of NGC 253 and NGC 4945 show ^-band at- 
tenuations of 30-100 magnitudes (Genzel et al. 1998). In 
a number of well-studied dwarf and irregular galaxies, a 
few compact regions dominate the infrared flux of the 
galaxy but remain invisible in the optical: Henizc 2-10 
(Kobulnicky & Johnson 1999; Vacca et al. 2002; Cabanac 
et al. 2005), NGC 5253 (Gorjian et al. 2001), SBS 0335- 
052 (Plante & Sauvage 2002; Dale et al. 2001; Houck et 
al. 2004), and II Zw 40 (Beck et al. 2002). Soifer et al. 
(2001) found that more than half of the infrared luminos- 
ity in infrared-luminous starburst galaxies is generated in 
regions 100 pc-1 kpc in size, and in the Antennae merger 
system, Mirabel et al. (1998) found that the most intense 
star-forming region is located outside both galaxy nuclei 
in an ~50 pc area that is entirely unremarkable at optical 

1* The attenuation at a given wavelength A is defined as Ax = 
k{X) X E{B — V) where fc(A) is the assumed reddening curve and 
E{B — V) is the color excess (e.g., Calzetti et al. 2001). For com- 
parison purposes, in this paper we assume Afjd = 0.82 X Ay- 
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wavelengths. 

In contrast to these highly obscured cases, most nor- 
mal H II regions show attenuations of only a few mag- 
nitudes. In the Large Magellanic Cloud, an Ay ~ 1 
magnitude is typical (Caplan & Deharveng 1986; Bell et 
al. 2002). Balmer decrement measurements for optically- 
selected H II regions in samples of spiral galaxies yield 
attenuations of Ay ~0.3-3.5 magnitudes (Zaritsky et al. 
1994). The mean F-band attenuation derived from ra- 
dio continuum measurements for H II regions in nearby 
galaxies is ~1. 7-2.0 magnitudes (Israel & Kennicutt 1980; 
van der Hulst et al. 1988). 

Probing to smaller scales in the Milky Way, we see ex- 
amples of a highly- embedded phase of star formation in 
the form of compact, ultracompact, or even "superultra- 
compact" H II regions which range in size from 1 pc down 
to 0.01 pc (Kurtz 2000). A study of the compact H ii re- 
gion S88B measured attenuations of Ay < 30 along with 
attenuation gradients that suggested gaps in the molec- 
ular cloud had allowed optical emission to emerge along- 
side the thermal dust emission in the infrared (Goetz et 
al. 2003). Results from carbon recombination line stud- 
ies of ultracompact H ii regions such as those in the W48 
star- forming complex are consistent with being pressure- 
confined (Roshi et al. 2005a, b). 

The new capabilities of the Spitzer Space Telescope 
(Werner et al. 2004) allow us to investigate the infrared 
properties of star formation within nearby galaxies down 
to scales as small as several hundred parsccs, the size of 
large H ii regions. We are thus able to ask the following 
question: what fraction of infrared-selected star-forming 
regions in normal galaxies are highly obscured? In other 
words, how much of the star formation in these galaxies 
have we been missing on these scales by only looking in 
the optical portion of the spectrum? 

Spitzer also gives us a new way to investigate radial at- 
tenuation properties of nearby galaxies. Previous stud- 
ies have found evidence for decreasing radial gradients 
in attenuation in samples of late-type galaxies using far- 
infrared to ultraviolet (FIR/UV) profiles and UV colors 
(Boissier et al. 2004, 2005; Popescu et al. 2005; Hill et al. 
1997) and in highly-inclined galaxies using broadband 
optical and near-infrared colors (Jansen et al. 1994). 
In a detailed study of NGC 300, Roussel et al. (2005) 
found a decreasing radial gradient in total infrared to 
far-ultraviolet (TIR/FUV). They suggest that while vari- 
ations in the metallicity, and presumably in the dust- 
to-gas ratio, may play a role in producing this gradient, 
other factors such as local variations in geometry or spec- 
tral energy distribution (SED) may be more important. 
In contrast, Zaritsky et al. (1994) found radial gradients 
in metallicity using optical spectra but saw only weak 
evidence for a gradient in attenuation. With the advent 
of Spitzer, these questions can now be investigated on 
the scale of individual star-forming regions. 

Taking advantage of the improved mid- and far- 
infrared resolution and sensitivity of Spitzer, the Spitzer 
Infrared Nearby Galaxies Survey (SINGS) is a multi- 
wavelength imaging and spectroscopic survey of 75 galax- 
ies within 30 Mpc (Kennicutt et al. 2003). Designed to 
span the range of properties seen in nearby, normal galax- 
ies, SINGS includes large spiral galaxies such as M51 
and M81 as well as dwarf galaxies. In addition to the 
Spitzer MIPS, IRAC, and IRS datasets, SINGS boasts 



a suite of ancillary data from radio continuum, HI, and 
CO to groimd-based BVRIJHK, Pa, Ha, UV, and X- 
ray wavelengths. This combination of high resolution 
and multiwavelength coverage opens up the possibility of 
studying the attenuation properties of infrared-selected 
star-forming regions drawn from the local population of 
nearby galaxies. 

In this work, we investigated the attenuation proper- 
ties of infrared-selected star-forming regions drawn from 
a sample of SINGS galaxies using multiband photome- 
try and an empirical method for deriving attenuations 
of infrared-selected star-forming regions. The paper is 
organized as follows. We describe our galaxy sample in 
Section 2 and the observation and reduction steps in Sec- 
tion 3. In Sections 4 and 5 we describe our method of 
multiband photometry and our treatment of Ha attenu- 
ation measurements. In Section 6 we present evidence of 
a paucity of large highly-obscured regions and of galaxy- 
wide radial trends in attenuation. We discuss the impli- 
cations of these results for the lifetimes of star-forming 
regions in Section 7, and we conclude in Section 8. 

2. THE GALAXY SAMPLE 

The SINGS sample was selected to span the range of 
properties seen in nearby normal galaxies. While not 
numerically representative of the local volume (Kenni- 
cutt et al. 2007b), the survey was constructed to contain 
objects with the full range of morphological types, lumi- 
nosities, and FIR/optical luminosity ratios (Kennicutt 
et al. 2003). The subset used in this work consists of 
38 galaxies and represents nearly all of the galaxies in 
SINGS that have discrete infrared star-forming regions. 
For reference. Table 1 gives the list of galaxies included in 
this work along with pertinent information. The quoted 
B magnitudes are total broadband magnitudes uncor- 
rected for Galactic extinction and inclination and derive 
from the Third Reference Catalog of Bright Galaxies (de 
Vaucouleurs et al. 1991, hereafter RC3). The galaxy di- 
mensions were taken from Kennicutt et al. (2003); the 
morphological types were taken from RC3. The inclina- 
tions were calculated using the quoted galaxy dimensions 
and position angles taken from the 2MASS Large Galaxy 
Atlas (Jarrett et al. 2003), and the following standard 
equation (Holmberg 1958): 

where a and b are the optical major and minor axis di- 
ameters, and where qo, the intrinsic optical axis ratio for 
an edge-on system, is set to 0.2 (TuUy & Fisher 1977). 
In the clearly edge-on cases - NGC 3190 and NGC 4631 
- the inclination was set to 90° , and the irregular galaxy 
NGC 6822 was set to an inclination of 0°. While inclina- 
tions derived in this manner suffer from systematic errors 
due to warps and spiral arms, kinematically-derived incli- 
nation measurements exist for only a subset of the galax- 
ies (The HI Nearby Galaxies Survey, THINGS: Walter 
et al. 2005). We have used the subset of our galaxy sam- 
ple for which kinematically-derived inclinations exist in 
order to estimate the error in applying equation 1. The 
quoted error is the standard deviation of three measure- 
ments: the inclination derived from 2MASS ellipticities, 
an inclination derived by fitting an ellipse to the IRAC 
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3.6 /im data, and the kincmatically-dcrivcd inclination 
from THINGS. Wc quote the average error (7 degrees) 
for all other cases in which kinematically-derived inclina- 
tions are not available. The distances were taken from a 
number of sources with preference given to direct dis- 
tance measurements from the literature; where direct 
measures were not available, we used flow-corrected esti- 
mates from Masters et al. (2007) (Hq = 70 km/s/Mpc). 

3. OBSERVATIONS AND REDUCTIONS 

The Spitzer images used in this work consist of the 
MIPS (Rieke et al. 2004) 24 fxm and the IRAQ (Fazio et 
al. 2004) 3.6, 4.5, 5.8, and 8.0 /xm datasets obtained as 
part of the SINGS project. The MIPS data reduction is 

described in Gordon et al. (2005), the IRAC data reduc- 
tion in Regan et al. (2004). The ground-based Ha and 
i?-band images were obtained at NOAO using the KPNO 
2.1-meter and CTIO 1.5-meter telescopes between 2001 
and 2003. A description of the SINGS observing strategy 
can be foimd in Kennicutt et al. (2003). 

Prior to beginning our analysis, the images were 
aligned, resampled to the same pixel scale, and con- 
volved to match the full-width-half-maximum (FWHM) 
of the MIPS 24 /im point-spread-function (PSF) using 
SINGS.TOOLi-\ Convolving to the 24 /im PSF miti- 
gates some of the inconsistency that would result from 
comparing higher resolution ground-based optical data 
with lower resolution 24 /im imaging where some frac- 
tion of the light from a given point source has spilled into 
neighboring regions. The 24 /im PSF, however, has fairly 
large diffraction lobes that are not replicated during our 
simple Gaussian convolution; for this reason, aperture 
corrections are applied at all wavelengths as described in 
Section 4. 

In order to correct the Ha images for the stellar 
continuum contribution, we subtracted off a scaled R- 
band image. The resultant images were corrected for 
[N ii]AA6548,6584A emission and the transmission curve 
of the Ha filter; these methods are described in more 
detail in Kennicutt et al. (2007b). Where available, the 
applied [N ii]/Ha ratios were taken from SINGS radial 
driftscan spectra, where available (Moustakas & Kenni- 
cutt 2006). Otherwise the corrections are approximate 
values derived from the B-band magnitude of the galaxy 
and the following empirical relation (Moustakas et al. 
2007): 

log{[N ii]A6584]/iJa) = -5.19 + 0A68Log{LB/LB ©) 

'(2) 

log([Ar ii]A6584]/ifa) = -0.187(Ms - 22.254) (3) 

where Mb,q = 5.47. This relation is based on a sample 
of integrated spectra of galaxies (Moustakas & Kennicutt 
2006; Jansen et al. 2000) with the AGN component re- 
moved as described in Moustakas et al. (2006). Table 2 
summarizes the correction factors used. 

4. REGION SELECTION AND MULTIBAND PHOTOMETRY 

Since most regions are unresolved at the 24 /xm reso- 
lution (~6 arcseconds), we selected regions in the 24 /xm 

^■^ SINGS-TOOL is an IDL program developed by Don Lindler 
(2006) for the SINGS Team for the automatic registration, align- 
ment, resampling, and PSF-matching of multiwavelength images. 



image using DAOfind, a point-source selection task in- 
cluded in the DAOphot stellar photometry package (Stet- 
son 1987). For the initial selection we set a detection 
threshold of lOcr, which corresponds to a typical infrared 
luminosity limit of ~3xl0^^ erg s~^); the infrared lumi- 
nosity limits for each galaxy are given in Table 1. For 
reference, the lowest corresponding Ha luminosity limit 
is comparable to the luminosity of a few Orion-type H ii 
regions (~ 3 x lO'^^erg s~^). The median 24 /xm lumi- 
nosity of the region sample is ~10^°, corresponding to a 
median Ha luminosity of ~2xl0^^ erg s~^ or a few tens 
of Orion-type H ii regions. 

We measured region fluxes using 500 pc diameter cir- 
cular apertures (~ 13" for a galaxy at 8 Mpc). While 
this aperture size is larger than a typical H II region 
(^100 pc), it was chosen as a compromise in order to ac- 
commodate the large 24 /im PSF and to allow a compar- 
ison of regions of the same physical size between galaxies 
at different distances. As an example, in Figure 1 we 
show the 24 fxm image of NGC 5194 with the final re- 
gion selection of 500 pc circular regions marked. In some 
cases, a 500 pc region was smaller than the FWHM of 
the 24 /im PSF; the aperture corrections discussed below 
were used to address this issue. 

Measuring robust fluxes for this work was complicated 
by the issue of how to determine the background value 
that should be subtracted from the raw region flux mea- 
surement. In addition to the sky background of the im- 
age, each 24 /im source sits on top of a diffuse component 
associated with the overall galaxy profile and internal 
structures such as spiral arms. Simple aperture photom- 
etry methods dictate using an annulus around each aper- 
ture to get a measure of the local background; however 
at the 24 /im resolution many star-forming regions lie 
within what would make up the sky annulus of a neigh- 
boring region, thereby contaminating the background es- 
timate. On the other hand, ignoring the contribution of 
the diffuse component and subtracting only a global sky 
background is not entirely satisfactory either as it sys- 
tematically overestimates the flux by 1-10% for regions 
on spiral arms or near the center of the galaxy. 

For this study, therefore, we employed three separate 
methods of sky subtraction. First, we used a global sky 
subtraction method by subtracting a simple sky back- 
ground defined as the a-clipped mean value from two 
off-galaxy portions of the image. Second, we employed 
a local background subtraction method where we mea- 
sured the background within large annuli (spanning from 
the measurement aperture to an outer radius equal to 4 
times the measurement aperture radius) centered on each 
region using a cr-clipped mode or mean estimation tech- 
nique. The large annuli are needed to minimize contam- 
ination due to the inclusion of neighboring IR-selected 
regions. To further reduce contamination, we removed 
pixels assigned to neighboring IR-selected regions prior 
to calculating the local background value. The success 
of this technique depends on the size of the background 
annuli used as well as properties of the galaxy such as 
the inclination, the spacing of the spiral arms, and the 
surface density of star-forming regions across the face of 
the galaxy. Figure 1 shows a small subset of the back- 
ground annuli used in the case of NGC 5194. Finally, we 
applied a profile background subtraction method which 
involved dividing the galaxy into elliptical annuli accord- 
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ing to the inclination and position angle and deriving a 
(7-clipped mean background value at each radius. As in 
the local background subtraction method, we removed 
pixels belonging to IR-selected regions prior to calculat- 
ing the background value in order to reduce the effects of 
contamination. This method is less sensitive to neighbor 
contamination than the local background method, but it 
does assume azimuthal symmetry. Fortunately, a major- 
ity of the targets in our current study are spiral galaxies 
that do not deviate far from this assumption. 

As an independent check on our background meth- 
ods, wc compared our flux measurements with those of 
Calzetti et al. (2005), a detailed multiwavelength study 
of NGC 5194. In that study background subtraction of 
region flux measurements was done using many rectangu- 
lar background areas defined by hand with statistics com- 
puted on the portion of the background area appropriate 
to the region being measured. Wc make the assumption 
that the level of fine-tuning demonstrated in Calzetti et 
al. (2005) for NGC 5194 provides the closest approxima- 
tion of the true local background for each region that 
one could achieve in treating each galaxy by hand. In 
the first three panels of Figure 2 we show a comparison 
between the 24 /um and Ha luminosities measured with 
each of the three background methods versus those of 
Calzetti et al. (2005) for a matched sample of 13 arsec- 
ond (-517 pc) regions in NGC 5194. We find that the 
most luminous regions are measured consistently regard- 
less of the background subtraction method because the 
background flux makes up a very small fraction of the 
raw region measurement. For faint, low surface bright- 
ness regions, however, the background is a more con- 
siderable portion of the measurement and thus the final 
value is highly dependent on the method applied. The 
global sky subtraction systematically underestimates the 
background because it neglects the diffuse component of 
the galaxy. The local background subtraction using large 
sky apertures is robust in many cases but tends to fail in 
crowded areas. The profile background subtraction tends 
to overestimate the background in interarm regions and 
underestimate the background in the spiral arms. Taken 
together, there is no standard method for dealing with 
the complexity of these nearby galaxies that is not in 
some way biased by the assumptions made or by the 
manner in which it is applied. In this work, we have 
taken the average of these three methods as our best es- 
timate of the flux and quote a standard error of the mean. 
In doing so, we reject individual measurements that are 
negative (i.e., the background was estimated to be higher 
than the source) as unphysical and average the remaining 
two measurements; this occured in only a small fraction 
(~5%) of the region sample. The success of the average 
versus the results of Calzetti ct al. (2005) is shown in the 
last panel of Figure 2. 

We estimated aperture corrections using the available 
MIPS and IRAC PSFs^^ as well as Gaussian PSF mod- 
els derived from measurements of the Ha and R images. 
Since the images are all convolved with a large Gaussian 
kernel in order to match the 24 jim data, the result- 
ing PSF profiles are nearly Gaussian, particularly in the 
ground-based data. However, the aperture corrections 

^® http://ssc.spitzer.caltech.edu/mips/psf.html; 
http: / / ssc.spitzer.caltech.edu/irac/psf.html 



for the ground-based and IRAC images are higher than 
they would be for unconvolved data due to the fact that 
our convolution does not precisely replicate the 24 ij,m 
diffraction lobes. The aperture corrections were indepen- 
dently derived for each galaxy depending on the aperture 
sizes used; Table 3 summarizes these corrections. 

Contaminant sources that show up in the 24 fiia images 
arise from either the Rayleigh- Jeans emission of fore- 
ground stars or the redshifted emission of background 
galaxies. We flagged objects that are located at a dis- 
tance greater than R25 from the center of the galaxy 
or that are included in the USNO catalog (Monet et 
al. 1998) or 2MASS point source catalog (Skrutskie et 
al. 2006). We also removed the 24 /im detections corre- 
sponding to the galaxy nucleus to avoid contamination of 
our region sample by embedded AGN. To verify the con- 
taminant selection and ensure that the 2MASS catalog 
did not accidentally identify objects within the galaxy, 
we then did a visual inspection of the images to confirm 
or reject the flagged objects as contaminants based on 
morphology in the unconvolved IRAC and ground-based 
images. 

Depending on the distance of the galaxy, and therefore 
the luminosity sensitivity limit, the output of DAOFIND 
may include many regions that, when measured with 
500 pc apertures, are signiflcantly overlapping. The 
net effect is that the same patch of a galaxy coiild be 
measured multiple times, and while this double-counting 
would not affect our results qualitatively, it would arti- 
ficially boost the number of 500 pc regions measured in 
the nearest galaxies. To limit this problem, we culled 
the region sample, weighting by a region's 24 /im flux in 
order to remove cases that overlapped spatially with a 
nearby, brighter region by more than 20%. 

5. Ha ATTENUATIONS 

Most previous measurements of dust attenuations for 
H II regions in normal galaxies have been based on 
Balmer decrements (e.g., Ha/H/3 ratios) derived from 
optical spectra (e.g., McCall et al. 1985; Caplan & De- 
harveng 1986; Zaritsky et al. 1994), or in a few instances 
Pa/Ha decrements (Quillen & Yukita 2001; Scoville et 
al. 2001). However such studies are based on star- 
forming regions selected from optical or near-infrared 
emission lines, and since our primary interest is in very 
dusty regions, we need a method for estimating the neb- 
ular attenuations that is valid for regions with up to 
several magnitudes of visual attenuation. A few au- 
thors have used multi-frequency radio continuum maps 
together with Ha photometry to estimate thermal radio 
fluxes and Ha attenuations, but this method is limited 
to the brightest 10 30 H 11 regions in the nearest galax- 
ies (e.g., Kaufman et al. 1987; van der Hulst et al. 1988). 
Instead we use the combination of mid-infrared luminosi- 
ties and Ha luminosities of the star-forming regions to 
derive attenuation estimates for an infrared-selected sam- 
ple. 

Our method is based on the recent work by Kennicutt 
et al. (2007a), who have generated a prescription for es- 
timating Ha attenuations {Ahu ) of star- forming regions 
using measurements of the 24 /xm and Ha fluxes for in- 
dividual star-forming regions. This empirically-derived 
method treats the 24 /xm flux of a region as a tracer 
of the obscured portion of the emission due to star for- 
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mation and the uncorrected Ha flux as a tracer of the 
unobscured portion. The attenuation is then related to 
the ratio of the total [24/xm+Ha] emission relative to the 
unobscured emission: 

Ahc. = 2.5 log[l + 0.038 vL^{2A)/L{Ha)], (4) 

where L(Ha) stands for the uncorrected Ha luminosity 
in erg s~^ and vL^(2'i) indicates the "monochromatic" 
24 /Ltm luminosity also in erg s~^. The coefficient 0.038 
in eq. 4 was calibrated using Ha/Pa attenuation mea- 
surements of 42 H II regions in M51 from Calzetti et al. 

(2005) . 

As is emphasized in Kennicutt ct al. (2007a), this 
method provides only an approximate estimate of the 
attenuation on an object by object basis, because the 
actual ratio of 24 /xm luminosity to Ha luminosity of 
a given object will depend on variable factors such as 
the geometry of the stars, dust, and gas and the ratio 
of ionizing flux to bolomctric luminosity of the OB as- 
sociation, which in turn varies with the age and mass 
function of the embedded cluster. This scatter can be 
quantified by comparing attenuations derived from eq. 4 
with independent measurements from the Pa/Ha decre- 
ments or other methods. Such comparisons have been 
carried out for large samples of H ii regions in M51 and 
M81 by Calzetti et al. (2005) and Perez- Gonzalez et al. 

(2006) and for 220 regions in 33 galaxies by Calzetti et al. 

(2007) . Comparisons of attenuation-corrected Ha fluxes 
with those derived via this infrared method show dis- 
persions of ±0.3 dex, which corresponds to uncertain- 
ties in the individual Ha attenuation measurements of 
±0.6—0.8 magnitudes. This confirms that this method 
provides only a crude estimate of the Ha attenuation 
for any individual H II region. However since the main 
aim of this paper is to identify strongly obscured regions 
{■^Ha ^ 3 magnitudes), this precision is more than ade- 
quate for our purposes. Nevertheless readers should bear 
in mind the uncertainties associated with individual at- 
tenuation measurements. 

6. RESULTS 

6.1. Distribution of Star-Forming Region Attenuations 

In order to address the question of how many infrared- 
selected regions are highly obscured, we used 24 /zm and 
Ha photometry to investigate the statistics of obscured 
star formation on 500 pc scales in a subset of SINGS 
galaxies. For our primary region sample, we identi- 
fied a total of 1772 regions using the 24 fMn. datasets 
as discussed above. In Figure 3a we show a histogram 
of Log[i/iy,/(24)/L(Ha)] along with the corresponding 
Ha attenuations. From the distribution of attenua- 
tions it is immediately clear that we do not see a large 
sub-population of highly-obscured star-forming regions 
within this infrared-selected sample on 500 pc spatial 
scales. The median value is 1.8 which corresponds to an 
Ha attenuation of 1.4 magnitudes (Ay ~ 1.7), and there 
is a slight tail to attenuations higher than 3-4 magni- 
tudes (skew = 0.35). Contaminant sources such as fore- 
ground stars or background galaxies would also popu- 
late the high attenuation portion of parameter space; to 
address this concern, we carefully verified that contam- 
inant sources had been removed from the sample Tising 
the approach described in Section 4. At this point, we 
must make clear what is meant by the "highly obscured 



population" or "subpopulation." Our beam encompasses 
many tens or hundreds of massive stars (see Table 1, 
minimum luminosities) and entire giant molecular clouds 
(GMCs). If all the stars are uniformly attenuated, then 
Ah a in Figure 3 represents the attenuation to each star. 
However, if, as is likely, some of the stars lie outside the 
CMC (due to breakout, blister Hii regions, etc.) and 
some of the stars are embedded, then the Ana repre- 
sents an average attenuation, weighted heavily toward 
the unattenuated stars. In this study, we do not find ev- 
idence for any large subpopulation of uniformly highly- 
obscured star- forming regions on 500 pc scales; we dis- 
cuss the question of smaller scales in Section 7.1.2. 

Figure 3b presents the same data in the form of a cu- 
mulative histogram. For illustration, the cumulative dis- 
tribution can be fit by a simple power law with an index 
of ~1.4 with high significance [R^ = 0.98) suggesting 
that the higher attenuation regions are simply an exten- 
sion of a smooth distribution. To quantify the fraction 
of regions in the highly-obscured tail relative to the en- 
tire sample, we make a cut at Log[i^Ly(24)/L(Ha)]= 2.7 
(A/fa ^ 3.3, A\/ ^ 4.0). This cut corresponds to ^2 
times the standard deviation of a Gaussian fit to the 
overall region distribution. By this definition, the num- 
ber of regions in the sample found in this highly-obscured 
tail of the distribution is 69, ~4% of the sample. 

Thus far, we have treated the entire sample of regions 
as one; however, the attenuation distribution changes 
slightly if we divide the sample l)y iiic;lination and lu- 
minosity. In Figure 4, we have split the sample into 
two panels by the inclination of the host galaxy, where 
high inclination is defined here as i > 60°. The median 
Log[j/L^(24)/L(Ha)] values from the low and high incli- 
nation subsets arc essentially equal within the errors — 
1.8 for the low inclination and 1.9 for the high inclina- 
tion cases. More striking, however, is the difference in 
shape between the two distributions: the low inclination 
subsample distribution is more sharply peaked (kurto- 
sis=0.71, with signifying a normal distribution) than 
that of the high inclination subsample (kurtosis=-0.28). 
This likely reflects the fact that the use of our empir- 
ical attenuation method becomes more complicated at 
high inclination; measured attenuations may be boosted 
in some cases due to increased pathlcngth through the 
disk, causing larger scatter. In addition, at the highest 
inclinations the region selection will be affected by con- 
fusion, and there may be a bias towards selecting more 
regions from the near side of the galaxy. In each panel of 
Figure 4, we divide the subsample into two bins of total 
galaxy 24 /xm luminosity, divided at the region-weighted 
mean galaxy 24 /xm luminosity of the sample (8.4x10^^ 
erg s"""^). The region attenuation distributions are sim- 
ilar to that of the full sample, with minimal change in 
the median attenuation as a function of galaxy luminos- 
ity. Thus, highly-obscured regions are drawn from both 
high and low luminosity galaxies. 

In this study we are interested in the attenuation prop- 
erties of star-forming regions on large, 500 pc scales. It 
is important to note, however, that this aperture may be 
much larger than the size of a given star-forming region 
as viewed at 24 /xm. To get a sense of how the chosen 
aperture size affects our results, we have done photome- 
try at scales spanning between roughly 200 and 1000 pc 
for most of the sample, between 100 and 1000 pc for 
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NGC 2403 (a particularly nearby case at 3.2 Mpc), and 
down to 20 pc for NGC 6822 (a dwarf galaxy at 0.5 Mpc). 
In Figure 5 we plot the median galaxy 24/zni/Ha ratio 
versus aperture size for a representative set of 4 galax- 
ies: NGC 0925, NGC 2403, NGC 4725, and NGC 6822. 
The vertical dashed line represents the point at which the 
aperture diameter is equal to the FWHM of the 24 /xm 
PSF, and below which the aperture measurement effec- 
tively undersamples the PSF. We see very little variation 
in the median attenuation down to 100 pc scales, while 
the scatter increases slightly over this range with decreas- 
ing aperture size. In our most nearby case, NGC 6822, 
we do see a significant rise in our measured attenuation 
as we probe below 100 pc, but at these scales we are also 
likely seeing the effects of the complex morphologies of 
Ha structures - shells and bubbles - that may be offset 
spatially from the brightest corresponding 24 /im soiucc. 
For the purposes of the current study, this analysis il- 
lustrates the fact that our application of the empirical 
method applies to large scales, i.e. greater than 100 pc, 
where the results are fairly insensitive to the particular 
aperture diameter chosen. 

6.2. Infrared Colors of Star-Forming Regions 

The infrared colors of the entire sample of regions are 
presented as flux density ratios in Figure 6. Follow- 
ing the discussion of Dopita ct al. (2005) we have cho- 
sen F(8.0)/F(4.5) versus F(24)/F(4.5) as a representative 
pair of colors that include contributions from the stellar 
(4.5 /im), PAH (8.0 //m), and dust (24 /;,m) components 
of a given region, although wc note that the IRAC bands 
do not sample the 10 /im silicate feature. have cor- 
rected the 8.0 /Ltm flux densities for the contribution of 
stellar emission using the following relation: 

F{8.0)corr = F{8.0) - 0.232 x F(3.6) (5) 

where F(8. 0)corr if' the corrected 8.0 /im flux density, 
F(8.0) is the observed value, and F(3.6) is the 3.6 /um 
stellar flux density (Helou et al. 2004). The contribution 
of dust emission to the 3.6 and 4.5 /im bands is assumed 
to be negligible for the purposes of this work. We have 
not applied any additional correction for the complex 
scattering that occurs when light from extended sources 
is incident on the IRAC focal plane^^ and that cause flux 
measurements to be systematically overestimated. The 
effect is negligible for the more distant galaxies in the 
sample, but even for the nearby cases the magnitude of 
the shift is less than 0.11 dex in Log[F(8.0)/F(4.5)] and 
less than 0.03 dex in Log[F(24)/F(4.5)]; this offset does 
not affect our conclusions. The contour encloses 95% of 
the total sample of regions, and the open and filled circles 
represent the individual regions and the highly-obscured 
subset, respectively. We can see that the highly-obscured 
regions do not stand out from the rest of the region sam- 
ple in terms of their infrared colors. The positive corre- 
lation seen between the F(8.0)/F(4.5) and F(24)/F(4.5) 
appears to plateau at Log[F(8.0)/F(4.5)]--1.25. This be- 
havior is seen in the region sample as a whole as well as 
in some of the individual galaxies. One plausible expla- 
nation is that at low fluxes, both the 8 /im PAH emission 
and the 24 /xm emission are excited by single optical/UV 
photons and are therefore proportional to the number 

^"^ http://ssc.spitzer.caltech.edu/irac/calib/extcal/ 



of absorbed stellar photons. This would imply that at 
low values of i^(24) one would expect a linear correlation 
of F(8.0) to i^(24), as is observed. As the 24 /xm emis- 
sion increases, however, larger grains could be sufficiently 
heated to provide a contribution to 24 fj,m emission from 
the Wien tail of this steady state warm dust. This then, 
may provide the plateau (increasing 24 yum emission ver- 
sus constant 8 /xm emission). 

In Figures 7-10, we plot the measured 24/im/Ha atten- 
uation versus the region 24 /xm luminosity, again using 
a 95% contour to represent the entire region sample; the 
dotted line represents the subsamplc defined by atten- 
uation in the previous section {Ana > 3.3, Ay > 4.0). 
In the entire sample, we see a general trend for more 
attenuated regions to be more luminous, a trend that 
is roughly reproduced within individual galaxies, with a 
few exceptions like NGC 3031. 

6.3. Properties of the Most Obscured Regions 

In applying the arbitrary cut from Section 6.1 [Ana > 
3.3, Av > 4.0), we were aiming to illustrate the relative 

paucity of uniformly highly obscured regions in the sam- 
ple (~4%), but we did not mean to imply that they are 
physicaUy distinct from the rest of the population. To 
demonstrate this fact we have further characterized the 
regions that make up the tail of the attenuation distri- 
bution (see Figure 3). In Table 1, wc list for each galaxy 
the total number of regions selected, the 24 /xm lumi- 
nosity of the faintest region measured, and the number 
of highly-obscured regions using our arbitrary definition. 
The number of such regions in a galaxy is typically small, 
and there are many galaxies with none. The exceptions 
are NGC 5713 and NGC 3198, in which all but a hand- 
ful of the measured regions qualify as highly obscured by 
our definition. As discussed in the next section, these 
two cases show globally high attenuation relative to the 
rest of the galaxy sample. 

The first question we ask is where the more highly- 
obscured regions are located within the host galaxy. In 
Figure 11 we present a histogram of the deprojected 
radius of the more highly-obscured regions within the 
host galaxy compared with the entire region sample (ex- 
cluding the edge-on cases NGC 4631 and NGC 3190). 
The histogram for the full sample extends out to 1.38 
R25 and has a median radius around 0.38 R25. In con- 
trast, the population in the highly-obscured tail is only 
found within 0.59 R25 with a median position of 0.26 
R25. Thus, the highly-obscured subset is drawn primar- 
ily from the inner parts of the host galaxies. To charac- 
terize these regions further, we investigate the colors of 
the highly-obscured tail compared with those of the rest 
of the regions in the host galaxy and with those of the 
entire region sample. Returning to Figure 6, we see that 
for the most part, the highly-obscured cases do not stand 
out in terms of infrared colors, but instead are scattered 
throughout the rest of the region color distribution of 
the host galaxy. In FiguK^s 7-10, w(^ s(;c that the loose 
trend between region luminosity and attenuation means 
that the regions which make up the tail of the attenua- 
tion distribution tend to also be more luminous. These 
comparisons all suggest that the more highly obscured 
regions in our sample are not a distinct sub-population; 
rather, they are simply the bright, central regions of their 
host galaxies. 
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6.4. Radial Attenuation Gradients 

Looking at the global properties of our galaxy sample, 

we find that a majority of the galaxies show a distinct 
and nearly linear decreasing radial trend in 24/im/Ha, 
a finding consistent with other studies (Boissicr ct al. 
2004, 2005; Popescu et al. 2005; Roussel et al. 2005). Fig- 
ures 12-16 illustrate this result by showing the 24/um/Ha 
ratio versus deprojected radius for each galaxy, with the 
galaxy sample split into two groups by inclination. The 
error bars are the standard errors of the mean derived 
using the three different background measurements de- 
scribed in Section 4. Galaxies that do not contain enough 
regions to determine the extent of a radial trend are not 
shown (NGC 1291, NGC 1316, NGC 1705, NGC 4450, 
Tololo 89, IC 4710); NGC 6822 is similarly excluded 
due to complicated morphology. The edge-on nature of 
NGC 4631 and NGC 3190 preclude an accurate radial 
deprojection. In cases where the region profile is well- 
described by linear regression we have shown the best 
linear fit to the radial range 0.0-0.65 i?25- In a few cases, 
the region profile is not adequately described by a linear 
fit due to more complex structure (e.g. bends in the 
profile); we will discuss this further in Section 7.2. 

In using the empirical conversion outlined in Section 
5, we are making the assumption that attenuation is the 
only mechanism responsible for varying the 24/Ltm/Ha 
ratio. We will proceed with this assmuption here but 
will return to the question of what else could be caus- 
ing the observed radial gradients in Section 7.2. Using 
an unweighted, linear least-squares fit to the radial pro- 
files we have investigated the relationship between these 
gradients and global galaxy properties of luminosity and 
inclination. We use the zeropoint, the linear extrapola- 
tion of the attenuation at the galaxy center, as a typical 
attenuation measure for each galaxy, and in Figure 17 we 
plot the best fit attenuation slope and zeropoint against 
galaxy Jf-band luminosity (Jarrett et al. 2003) and in- 
clination. We find no strong trends; a tentative correla- 
tion appears to exist between the attenuation zeropoint 
and the galaxy H-hand luminosity and is corroborated 
by a similar trend between Mb and zeropoint. Dale et 
al. (2007) studied the global total infrared-to-ultraviolet 
(TIR/UV) ratios, often used as a proxy for attenuation, 
within the entire SINGS sample and found that TIR/UV 
increases with increasing iJ-band luminosity. Our data 
are consistent with this result, though the relationship 
between attenuation zeropoint and galaxy iJ-band lu- 
minosity is not as strong. The galaxy inclination does 
not appear to regulate either the observed slope or ze- 
ropoint of the attenuation profile, consistent with find- 
ings by Dale et al. (2007) on global TIR/UV ratios. In 
terms of morphology, we see a weak tendency for early 
type galaxies with bars to show somewhat steeper slopes 
that early type galaxies without bars; no such relation 
is seen for late typ(«. We note that most of the highly- 
obscured regions in our sample are found in barred galax- 
ies, consistent with other studies that have found an as- 
sociation between bars and circumnuclear star-forming 
regions (Regan et al. 2006; Kennicutt 1998). 

In the two overlapping cases between this work and 
Boissier et al. (2004) (NGC 5194 and NGC 3031) we 
find similar results. Boissier et al. (2004) used FUV/IR 
ratios to derive UV attenuation estimates and compared 



these with Ana measured from Balmer decrements of 
H II regions. In NGC 5194, they see a radial gradient of 
roughly 1 magnitude in Ha attenuation over the optical 
radius of the galax;y, compared with -^0.9 magnitudes in 
our study. In NGC 3031, they see a much shallower at- 
tenuation profile with a radial drop of ~0.1 magnitudes, 
comparable to our result (A^f^ ^0.2 magnitudes). 

Zaritsky et al. (1994) looked at the subject of attenu- 
ation gradients using optical spectra of a sample of 159 
H II regions from 14 galaxies. They found only weak 
evidence for a radial trend in attenuation (see their Fig- 
ure 7c). However, we can interpret this result within the 
context of our study by noting that due to their more 
limited sample, they plotted the Ay of the entire region 
sample versus fractional isophotal radius. Since we find 
that the attenuation zeropoint can vary from galaxy to 
galaxy, viewing the region sample as a whole tends to 
wash out a radial signature. Figure 18a shows our version 
of their Figure 7c; a radial gradient is evident, although 
it is much less clear than what is seen in any individual 
galaxy. When we randomly select 159 regions from our 
sample to match the sample size of Zaritsky et al. (1994), 
we obtain a result that is qualitatively indistinguishable 
from Zaritsky et al. (1994) Figure 7c (Figure 18b). 

In view of the range in attenuation zeropoints found 
within the sample and the fact that our targets display 
decreasing radial attenuation gradients, we can return to 
the question of the nature of the tail of the attenuation 
distribution. In defining a global attenuation cut, we 
selected the central regions of galaxies with the highest 
attenuation zeropoints. To visualize that this is indeed 
the case, we have included in Figures 12-16 a dotted line 
corresponding to the attenuation cut applied in Section 
6.1. Regions that fall above the dotted line arc prefer- 
entially the central regions of galaxies with overall high 
attenuation. This discussion underscores the fact that 
the highly-obscured fraction we quoted initially is quite 
arbitrary and dependent on the particular galaxy sam- 
ple. Nonetheless, while the precise value of the highly- 
obscured fraction is dependent on the sample, the state- 
ment that only ~4% of the regions in this particular set 
of galaxies have 24/im/HQ; attenuations higher than 2.7 
{Aho ^ 3.3) is representative of the robust result that 
these galaxies are not hosting a large number of 500 pc 
regions in which the vast majority of luminous stars are 
highly obscured. 

In a few galaxies there are a handful of regions that 
stand out relative to the rest of the regions in the host 
galaxy. These regions arc not necessarily liighly-obsc;urcd 
in a global sense - that is, they are not part of the highly- 
obscured tail of the global attenuation distribution - but 
they do have higher attenuations than the rest of the 
region population of the host galaxy near that radius. 
NGC 3938 and NGC 5474 arc examples, with a few of 
these local outlier regions at 0.6-0.8 R25. Future work 
will be needed to look into the question of whether these 
regions are in some way exceptional or whether they sim- 
ply reflect the stochasticity of the disk at large radii. 

7. DISCUSSION 

7.1. Obscured Star Forming Regions 

The main purpose of this study was to characterize the 
attenuation distribution of infrared-selected star-forming 
regions as measured on large scales. In other words, the 
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goal was to determine what fraction of bright infrared 
sources have Ha counterparts and what fraction could 
be defined as highly obscured. The median attenuation 
we found {Ana ~l-4 magnitudes) indicates that roughly 
72% of the Ha emission is typically "missing" from a 
500 pc region, but cases where the Ha is attenuated by 
more than 3 magnitudes are quite rare (^4% of our sam- 
ple). A complementary picture has been seen in direct 
measurements of mid-infrared silicate absorption within 
kiloparsec-sized regions in normal galaxies where most 
do not show foreground dust screen attenuation greater 
than Av ^ 3 (Smith et al. 2006). Similarly, optical nu- 
clear spectra of SINGS galaxies show an average Ay ^ 1 
(Dale et al. 2006). However, the fact that we see such a 
small fraction of uniformly highly-obscured regions could 
be due to a number of possible scenarios. 

7.1.1. Timescales 

While dramatic examples of highly-obscured sources 
do exist in the local Universe on a variety of size scales, 
the main result of this study is that the fraction of such 

regions in normal, star-forming disk galaxies is small on 
500 pc scales. Equivalently, we can say that the vast 
majority of current star formation in these galaxies is 
not hidden from view at optical wavelengths. If we were 
to make the very simple assumptions that all large star- 
forming regions go through a highly-embedded phase and 
that a single star cluster dominates on 500 pc scales, 
we can estimate the lifetime of such an obscured period. 
Taking a typical H II region lifetime to be 5-10 Myr, our 
observed obscured fraction of ~4% implies that a region 
would spend only ^ 10''' years being obscured by more 
than 3.3 magnitudes in Ha (Ay ^ 4 magnitudes). One 
could also account for these statistics if only a small frac- 
tion of the star-forming regions undergo a more lengthy 
obscured phase. These short time scales are in rough 
agreement with those inferred for the lifetime and dy- 
namical timescales of Young Stellar Object (YSO) out- 
flows in nearby Galactic molecular clouds, (e.g., Quillen 
et al. 2005; Zhang et al. 2005; Bontemps et al. 1996). Fur- 
thermore, our sample consists of large 500 pc regions per- 
haps containing hundreds of stars which should shorten 
the time for sufficient cloud disruption even further. 

That said, in the above estimate we have envisioned 
simplified regions in which all stars are uniformly ob- 
scured. In reality, each 500 pc region potentially encom- 
passes many highly-obscured, compact infrared knots as 
well as emerged clusters which, when blended together in 
the large beam, would yield a lower average attenuation. 
We will return to this issue in the following section, but 
we note here that the observed highly-obscured fraction 
could also be explained by envisioning a toy model in 
which each 500 pc region half of the massive stars are 
embedded while half have broken out of one side of the 
giant molecular cloud in which they formed. Roughly 
half of the emerged cases have broken out on the oppo- 
site side of the cloud, and thus remain obscured from the 
point of view of an observer over the short lifetime of 
massive stars. The net Ha absorption in this case would 
be of order 75% (Ajj„ = 1.5), similar to the measured 
median attenuation in our sample. This suggests break- 
out times of order one half of the massive star lifetime, 
or roughly 1 Myr. In either scenario, the small fraction 
of regions with nearly all of their stars highly-obscured 



indicates that breakout times must be less than the typi- 
cal lifetime of a star-forming region and also less than or 
of order of the main sequence lifetime of the dominant 
sources of luminosity. 

Although we have shown that infrared-luminous star- 
forming regions with relatively faint optical counterparts 
are very rare in the extended disks of the SINGS spirals, 
we have also seen that they are considerably more com- 
mon in the circumnuclear regions of the galax;ies. This 
difference may reflect the higher interstellar pressures in 
the nuclear regions, which will act to slow the expan- 
sion of H II regions and prolong the lifetime of any dust- 
embedded phase, as modeled, for example, by Dopita et 
al. (2005). These regions are also characterized by much 
larger mean gas column densities, however, so we would 
expect higher dust opacities even without the dynamical 
effects of the higher pressure environment. More detailed 
studies of individual regions will be needed to separate 
these effects. 

The strong association that we find between 24 fim 
sources and optical H ii regions also provides interesting 
contraints on the maximum lifetimes of the star form- 
ing clouds. If the molecular/ atomic clouds and their as- 
sociated dust remained intact for lifetimes much longer 
than that of the H II regions, we would expect many 
of the 24/xm-bright H li regions to be associated with 
the older clusters and OB associations, coincident with 
ultraviolet-bright regions. Since the dust can be heated 
by non-ionizing stars such regions would stand out as 
having anomalously high 24/xm/Ha ratios, effectively 
masquerading as highly obscured regions. While the ob- 
scured fraction that we find (~4%) is therefore an upper 
limit, the scarcity of regions with high 24/xm/Ha ratios 
means that 24/xm-bright clusters older than 5-10 Myr 
are even more rare. This suggests that the dust and gas 
clouds surrounding the star clusters/associations must 
largely dissipate on this time scale. Independent evi- 
dence for such time scales comes from observations of 
expansion lifetimes of HI giant shells in the LMC and 
SMC, which also show an absence of objects with expan- 
sion ages of more than 10 Myr (e.g., Kim et al. 1999; 
Hatzidimitriou et al. 2005). 

7.1.2. Size scales 

Our results do not preclude the presence of embed- 
ded protostars or ultra-compact H II regions within the 

larger complexes that we have studied. In this work 
we have presented attenuation measurements for large, 
500 pc apertures. A common aperture size allowed us 
to compare the same spatial scale between galaxies at 
distances ranging from 3 to 25 Mpc; however, a 500 pc 
aperture corresponds to a larger scale than that of a typ- 
ical H II region, let alone a compact H ii region (~1 pc). 
In every aperture applied to the 24 /xm data, we are 
potentially summing together a number of smaller in- 
frared star-forming knots with varying attenuations; the 
net result of this blending is that the attenuation mea- 
sured is an average for the region. This situation is of 
course only exacerbated in global star formation studies 
at high redshift in which measurements are necessarily 
averaged over an entire galaxy. In Section 6.1, we ad- 
dressed the question of whether probing to smaller scales 
would have drastically changed the fraction of highly ob- 
scured sources. Based on these tests we do not expect 
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our results on the highly-obscured fraction to change dra- 
matically down to 100 pc scales. When we repeat our 
analysis on similar datasets for the Antennae galaxies, for 
example, we are able recover the highly-obscured case de- 
spite the blending that occurs due to our large aperture. 
The region itself is ~50 pc with an attenuation of ~9-10 
magnitudes in V derived from NIR spectra (Gilbert et 
al. 2000); our empirically-calibrated photometric atten- 
uation measurement is ~4.8 in V {Ahu ~ 3.9) using a 
500 pc aperture, with the region being a clear (8cr) out- 
lier from the rest of the regions in the system. Thus, if 
a large sub-population of similarly luminous and embed- 
ded regions existed in our sample of SINGS galaxies, we 
should have found them in this study. 

This is not to suggest the absence of much smaller or 
less luminous embedded infrared sources. We know from 
stiidics of the Milky Way that ultracompact and compact 
H II regions exist on scales of ^0.1 - 1 pc with TIR lumi- 
nosities of - lO-'^^-lO-'^'' erg (Crowther 2005; Akabane 
& Kuno 2005). Star-forming regions in our sample are 
typically at least ~100 times more luminous than a com- 
pact H II region just in the 24/im-band lO'"*^ — 10'*^ erg 
s~^). No doubt a multitude of such regions exist within 
each galaxy in our sample, but they are not large enough 
or luminous enough to be studied with the current data. 

7.2. Radial Gradients 

The fact that we see few highly-obscured cases and 
such smooth radial gradients suggests that the measured 
attenuations on these scales are dictated to a greater 
extent by the global properties of the galaxy such as 
metallicity, gas surface density, inclination, geometry, or 
gradual changes in the region SED than by region-to- 
region fluctuations. Since metallicity is related to dust 
content and since decreasing gradients in metallicity are 
observed in many galaxies, it is reasonable to expect it to 
contribute to the attenuation properties of star-forming 
regions in a manner that varies with radius (Zaritsky 
et al. 1994). In contrast to previous work (Boissicr et 
al. 2004, 2005), a preliminary comparison of the metal- 
licity and attenuation gradients for the low inclination 
galaxies in our sample did not suggest any correlation 
between the two. However, this is perhaps not surprising 
as our sample does not cover a wide baseline in metal- 
licity, and gradients in gas column density in galaxies 
may play a larger role in the global attenuation profiles. 
Roussel ct al. (2005) found a similar radial gradient in 
FUV/TIR in NGC 300 and noted that while the disper- 
sion in FUV/TIR ratios of individual regions was at least 
a factor of ~2 (>0.30 dex), the standard deviation of 
abundances was only 0.11 dex. From this they suggested 
that metallicity was likely not the only factor responsi- 
ble for the gradient. The inclination of the galaxy is a 
consideration since distributed dust seen along a line-of- 
sight through the disk may cause additional attenuation 
of the region's Ha flux. In this sample, however, galaxy 
inclination does not appear to play a dominant role in 
regulating either the slope or the zeropoint of the global 
attenuation profile. Instead, the situation appears to be 
more complicated and may depend strongly on properties 
of an individual galaxy such as the thickness and mor- 
phology of dust lanes and the dust and gas scale heights. 

It is interesting to note that while the majority of the 
sample show linear radial trends in 24/im/HQ!, there are 



two cases (NGC 1097 and NGC 4254) that show a bend 
in the profile at roughly 0.6 R25. Since metallicity pre- 
sumably plays a role in the dust-to-gas ratio, it is rea- 
sonable to ask whether a bend in the attenuation profile 
could be related to a bend in the abundance profile of 
a galaxy. Previous work by numerous authors presented 
examples of bends in abundance gradients and discussed 
the question of whether they were real or simply an ar- 
tifact of the method used to determine (0/H) (Zaritsky 
et al. 1994, and references therein). Pilyugin (2003) ar- 
gued strongly that the bends were in fact artificial, that 
they were due to either an inappropriate application to 
low metallicity regions of the high-metallicity relation- 
ship between oxygen abundace and strong lines or due to 
a gradient in the excitation parameter within the disk. 
While beyond the scope of this paper, it is interesting to 
note that we find examples of bends in attenuation gra- 
dients using a completely different technique. This leads 
to the question of whether these bends could be related 
to gradients in the [N Il]/Ha ratio or systematics inher- 
ent to the empirical method of estimating attenuation, or 
whether they represent a real change in the attenuation 
properties of the disk. 

The results summarized above assume that all of the 
radial trends observed in the 24/Ltm/Ha ratios of the H ii 
regions arise from changes in dust attenuation. However 
it is possible that some of the radial trends that we ob- 
serve are produced by systematic changes in the infrared 
SED of the emitting regions with galactocentric radius, 
independent of extinction. For example, if the diame- 
ter of H II regions becomes systematically larger with 
galactocentric radius, this would lower the characteristic 
intensity of the heating radiation field at larger radii and 
tend to lower the mean dust temperatures and the frac- 
tion of total dust luminosity emitted in the 24-/xm band. 
Spatially-resolved measurements of the infrared spectral 
energy distributions of SINGS galaxies show little evi- 
dence for such radial gradients in infrared SED shape 
(Bendo et al. 2006; Dale et al. 2005; Perez-Gonzalez et 
al. 2006), but we cannot exclude the possibility that there 
is at least some small contribution from this effect to the 
gradients that we observe in Figures 12-16. 

8. CONCLUSIONS 

In this work, we have studied the incidence of obscured 
star formation and the distribution of attenuation prop- 
erties for infrared-selected star-forming regions measured 
on large scales (~500 pc) in a sample of 38 SINGS galax- 
ies. We have made use of an empirical method for deriv- 
ing Ha attenuations which combines the 24 /xm and Ha 
datasets (Kennicutt et al. 2007a). Under the assump- 
tion that the 24/im/Ha ratio can be used as a proxy for 
attenuation, our main results are as follows: 

1. We do not find a distinct sub-population of uni- 
formly highly-obscured regions on large (~500 pc) 
scales. Instead our data suggest a continuous dis- 
tribution of attenuations (power law index of ~1.4) 
within the entire sample of 1772 regions. 

2. With a few exceptions, those regions that do have 
higher obscuration {Ahc > 3.3) are the bright, cen- 
tral regions of luminous galaxies with high global 
attenuation. The regions in the tail of the attenua- 
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tion distribution represent a small fraction (~4%) 
of the region population. 

3. Most galaxies show a decreasing radial trend in 
24/Ltm/Ha; the current sample shows a range of 
slopes and zeropoints. A weak correlation may 
exist between the attenuation zcropoint and the 
galaxy luminosity, consistent with previous stud- 
ies, but the galaxy inclination does not appear to 
play a dominant role in either the observed slope 
or zeropoint of the profile. 

4. This study reports on the attenuation (Affa) prop- 
erties of large (~500 pc) star-forming regions, but 
does not comment on the prevalence of embed- 
ded infrared sources that no doubt exist on smaller 
scales and at lower luminosities similar to compact 
and ultracompact Hll regions in the Milky Way. 

5. A more general conclusion of our work is the high 

degree of spatial correlation between 24 /im and 
Ha emission in SINGS galaxies. In agreement with 
other studies (e.g., Perez-Gonzalez et al. 2006), this 
work suggests that for every bright 24 /xm source in 
a normal galaxy there is a visible Ha counterpart, 
a fact which opens the door not only to attenuation 
studies but also to the possibility of using the em- 
pirical [24/^m+Ha] combination as a new measure 
of star-formation rates in galaxies. 
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TABLE 1 

Galaxy Sample 



Galaxy 


Distance 


Ms 


D25 


Type 


i 




LogKi,(24)] 




Slope 


Zcropoint 














Limit 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 



84± 7 


7 


39.17±38.52 









53± 7 


11 


40.48±39.94 





-1.21±0.47 


2.37±0.19 


26± 6 


70 


39.46±38.40 


2 


-0.53±0.20 


2.51±0.07 


58± 4 


35 


39.24±37.95 


1 


-0.42±0.25 


2.47±0.09 


49± 7 


66 


39.50±38.29 


3 


-2.25±0.22 


2.73±0.08 


35± 7 


1 


39.30±37.52 









46± 7 


1 


40.94±40.10 









52± 7 


9 


39.77±38.87 





-1.52±0.21 


2.57±0.05 


38± 7 


62 


39.61±39.17 





-0.98±0.20 


1.81±0.07 


44± 7 


2 


39.77±38.91 









58± 7 


88 


38.36±37.18 





-0.47±0.16 


1.45±0.05 


67± 1 


28 


39.53±38.87 





-0.44±0.33 


2.12±0.14 


65± 3 


14 


39.51±39.18 





-0.14±0.33 


1.56±0.13 


64± 4 


81 


38.64±37.32 





-0.28±0.29 


1.43±0.13 


22± 3 


80 


39.40±38.77 





-0.63±0.15 


2.18±0.06 


90± 7 


5 


40.57±39.71 









70± 3 


43 


39.50±38.28 


35 


-0.82±0.28 


3.22±0.10 


69± 3 


13 


38.59±37.43 









49± 4 


13 


39.46±38.79 





-1.25±0.32 


2.30±0.12 


65± 8 


100 


39.61±39.08 





-1.56±0.13 


2.41±0.04 


65± 4 


71 


39.57±39.26 


2 


0.11±0.23 


2.17±0.09 


25± 7 


77 


39.46±38.66 





-0.96±0.12 


2.15±0.04 


30± 7 


119 


39.69±39.19 





-1.13±0.10 


2.03±0.04 


32± 7 


87 


39.64±39.40 





-0.87±0.13 


2.36±0.05 


42± 7 


3 


39.89±39.03 









68± 7 


25 


39.70±38.84 


1 


-0.80±0.54 


2.45±0.19 


90± 7 


68 


39.59±39.00 


1 






46± 7 


32 


39.44±38.96 





-2.16±0.89 


2.58±0.36 


54±12 


217 


39.33±39.08 





-0.38±0.08 


1.98±0.03 


54± 7 


3 


39.96±38.83 









27± 7 


15 


38.90±37.78 





0.01±0.37 


0.97±0.16 


27± 7 


14 


40.68±40.41 


13 


-1.74±0.53 


3.52±0.16 


40± 7 


5 


39.44±38.37 









0± 7 


13 


37.70±37.45 









32± 9 


158 


39.51±38.94 


4 


-1.02±0.13 


2.38±0.05 


73± 2 


77 


39.86±38.78 


7 


-2.15±0.27 


2.83±0.09 


38± 7 


11 


40.47±40.05 





-1.55±0.70 


2.50±0.29 


49± 3 


48 


39.14±37.80 





-0.76±0.17 


1.55±0.07 



NGC 0024 
NGC 0337 
NGC 0628 
NGC 0925 
NGC 1097 
NGC 1291 
NGC 1316 
NGC 1512 
NGC 1566 
NGC 1705 
NGC 2403 
NGC 2841 
NGC 2976 
NGC 3031 
NGC 3184 
NGC 3190 
NGC 3198 

IC 2574 
NGC 3351 
NGC 3521 
NGC 3627 
NGC 3938 
NGC 4254 
NGC 4321 
NGC 4450 
NGC 4536 
NGC 4631 
NGC 4725 
NGC 5194 
Tololo 89 
NGC 5474 
NGC 5713 

IC 4710 
NGC 6822 
NGC 6946 
NGC 7331 
NGC 7552 
NGC 7793 



7.3±2.3 (i) 
22.4±2.3 (i) 
7.3±1.5 (k) 
9.1±0.2 (b) 
17.1±2.3 (i) 
10.8±2.3 (i) 
24.3±2.3 (i) 
11.8±2.3 (i) 
20.3±2.3 (i) 
5.1±0.6 (m) 
3.2±0.2 (b) 
14.1±1.5 (h) 
3.6±0.4 (d) 
3.6±0.1 (b) 
11.1±2.0 (g) 
20.9±2.3 (1) 
13.7±0.5 (b) 
4.0±0.4 (e) 
9.3±0.4 (b) 
10.1±2.3 (i) 
9.4±0.1 (b) 
13.3±2.3 (i) 
16.6±2.0 (i*) 
14.3±0.5 (1) 
16.6±2.0 (i*) 
14.4±0.3 (b) 
8.1±2.3 (i) 
11.9±0.3 (b) 
7.8±2.3 (i) 
16.7±2.3 (i) 
6.8±1.4 (a) 
29.4±2.3 (i) 
9.0±2.3 (1) 
0.5±0.0 (j) 
6.8±1.4 (c) 
14.5±0.6 (b) 
21.0±2.3 (i) 
3.8±0.4 (f) 



-17.14±0.13 
-19.69±0.13 
-19.37±0.10 
-19.11±0.11 
-20.93±0.07 
-20.78±0.04 
-22.51±0.08 
-19.23±0.10 
-21.21±0.03 
-15.77±0.13 
-18.61±0.07 
-20.66±0.10 
-16.94±0.13 
-19.91±0.03 
-19.87±0.10 
-19.48±0.10 
-19.81±0.10 
-17.22±0.19 
-19.32±0.10 
-20.18±0.10 
-20.21±0.13 
-19.73±0.10 
-20.66±0.08 
-20.73±0.08 
-20.20±0.08 
-19.64±0.08 
-19.78±0.16 
-20.27±0.13 
-20.49±0.06 
-18.33±0.17 
-17.88±0.15 
-20.50±0.15 
-17.26±0.20 
-14.03±0.06 
-19.55±0.10 
-20.46±0.10 
-20.36±0.13 
-18.28±0.05 



5.8x 1.3 
2.9x 1.8 
10.5X 9.5 
10.5X 5.9 
9.3x 6.3 
9.8x 8.1 
12.0X 8.5 
8.9x 5.6 
8.3x 6.6 

I. 9x 1.4 
21.9x12.3 

8.1X 3.5 
5.9x 2.7 
29.6x14.1 
7.4x 6.9 
4.4x 1.5 
8.5x 3.3 
13.2X 5.4 
7.4x 5.0 

II. Ox 5.1 
9.1X 4.2 
5.4x 4.9 
5.4x 4.7 
7.4x 6.3 
5.2x 3.9 
7.6x 3.2 
15.5X 2.7 
10.7X 7.6 
11.2X 6.9 
2.8x 1.7 
4.8x 4.3 
2.8x 2.5 
3.6x 2.8 
15.5x13.5 
11.5X 9.8 
10.5X 3.7 
3.4x 2.7 
9.3x 6.3 



SA{s)c 
SB(s)d 
SA(s)c 
SAB(s)d 
(R)SB(r)b 
{R)SB(l)0/a 
(R)SAB(s)0 
SB(r)ab 
(R)SAB(rs)bc 
SAO pec 
SAB(s)cd 
SA(r)b 
SAc pec 
SA(s)ab 
SAB(rs)cd 
SA(s)a pec 
SB(rs)c 
SAB{s)m 

SB(r)b 
SAB(rs)bc 
SAB(s)b 
SA(s)c 
SA(s)c 
SAB(s)bc 
SA(s)ab 
SAB{rs)bc 
SB{s)d 
SAB(r)ab pec 
SA(s)bc pec 

SB(rs)dm 
SA(s)cd pec 
SAB(rs)bc pec 
SB(s)m 
IB(s)m 
SAB(rs)cd 

SA(s)b 
(R)SB(s)ab 
SA(s)d 



References. — (a) Drozdovsky & Karachentsev (2000); (b) Freedman et al. (2001); (c) Karachentsev et al. (2000); (d) Karachentsev et al. (2002); 
(e) Karachentsev et al. (2003a); (f) Karachentsev et al. (2003b); (g) Leonard et al. (2002); (h) Macri et al. (2001); (i) Masters et al. (2007); (j) 
Pietrzynski et al. (2004); (k) Shapley et al. (2001); (1) Tonry et al. (2001); (m) Tosi et al. (2001) 

Note. — Col. (1) Galaxy name; Col. (2) Distance (Mpc) with associated reference; Col. (3) Absolute B magnitude from de Vaucouleurs et al. 
(1991) (RC3); Col. (4) Major and minor axis diameters at the 25 magnitude isophote from RC3 (arcmin); Col. (5) Morphological type from RC3; 
Col. (6) Inclination derived using eq. 1 (degrees); Col. (7) Number of infrared-selected regions included in this study; Col. (8) Log of the 24^m 
luminosity of the faintest region (erg s~^); Col. (9) Number of highly obscured regions, as defined in Section 6.1; Col (10) Slope of radial attenuation 
gradient as defined in Section 6.4 (dex/kpc); Col (11) Zeropoint of radial attenuation gradient as defined in Section 6.4 (dex). 

Distance set to Virgo Cluster center. 



Highly-Obscured Star-Forming Regions in SINGS 

TABLE 2 
Ha Correction Factors 



[N ll] A6584/Ha Final Correction 



Galaxy 


Ratio* 




Factor 




NGC 0024 


0.268±0.008 


0, 


,8947±0, 


,0028 


NGC 0337 


0.197±0.002 


0, 


,9819±0, 


,0002 


NGC 0628 


0.497±0.014 


0, 


,9645±0, 


,0009 


NGC 0925 


0.233±0.006 


0, 


,7814±0, 


,0035 


NGC 1097 


0.485±0.004 


0, 


,7323d=0. 


,0017 


NGC 1291 


1.060±0.157 


0. 


,9487±0. 


,0051 


NGC 1316 


2.489±0.126 


0. 


,6646±0. 


,0103 


NGC 1512 


0.481±0.012 


0, 


,9713±0. 


,0005 


NGC 1566 


0.507±0.004 


0, 


,9737±0, 


,0002 


NGC 1705 


0.065±0.002 


0, 


,9560±0, 


,0012 


NGC 2403 


0.231±0.002 


0, 


.7679±0, 


,0012 


NGC 2841 


0.866±0.027 


0. 


.4972±0. 


,0059 


NGC 2976 


0.269±0.004 


0. 


.7424±0. 


,0020 


NGC 3031 


0.36±0.30'' 




0.67±0. 


,12 


NGC 3184 


0.439±0.009 


0. 


.6571±0. 


,0036 


NGC 3190 


0.958±0.091 


0. 


.4594±0. 


,0206 


NGU 3198 


U.34o±0.008 


0. 


.7125±0. 


,0036 


iC 2574 


U.ll±U.oU 




0.87±0. 


,18 


■\TOl^ QQC1 
INVjrL' OOOi 


n yicQ_L.r\ t\c\A 


0. 


.Oz9z±U. 


,0017 




u . 4Z D in u . uu'i: 


u. 




,UU1D 


NGC 3627 




0, 


6381 +0 


,0018 


NGC 3938 


0.310±0.009 


0, 


,7460±0, 


,0040 


NGC 4254 


0.307±0.002 


0, 


,7096±0, 


,0012 


NGC 4321 


0.396±0.005 


0, 


,8285±0, 


,0014 


NGC 4450 


0.873±0.043 


0, 


,4777±0, 


,0095 


NGC 4536 


0.432±0.005 


0, 


,8060±0, 


,0015 


NGC 4631 


0.34±0.30'' 




0.71±0, 


,12 


NGC 4725 


1.261±0.077 


0, 


,5627±0, 


,0105 


NGC 5194 


0.524±0.004 


0, 


,6023±0, 


,0013 


Tololo 89 


0.18±0.30'' 




0.90±0, 


,09 


NGC 5474 


0.15±0.30'' 




0.83±0, 


,17 


NGC 5713 


0.418±0.003 


0, 


,6565±0, 


,0011 


IC 4710 


0.12±0.30'' 




0.99±0, 


,01 


NGC 6822 


0.03±0.30t' 




l.OOiO, 


,03 


NGC 6946 


O.SliO.SOt" 




0.71±0, 


,13 


NGC 7331 


0.507±0.007 


0, 


,6433±0, 


,0025 


NGC 7552 


0.44±0.30'' 




0.97±0, 


,01 


NGC 7793 


0.287±0.005 


0, 


.9822±0, 


,0002 



Note. — The final correction factor multiplies 
the raw flux measurement and corrects for the 
[N ii] AA(6584,6548) /Ha ratios as well as the Ha filter 
transmission. Quoted errors are formal measurement 
errors only and do not include physical uncertainties, 
e.g., due to variations in the [N ii] AA(6584,6548)/Ha 
ratio within the galaxy. 

^[N ii] AA(6584,6548)/Ha ratios were de- 
rived as discussed in Section 3. The assumed 
[N ii] A6548/[N ii] A6584 ratio is 0.335, consistent 
with other studies (Calzctti ct al. 2005). 

^Ratio derived using cq. 2. The standard deviation 
of this relation is estimated to be ±0.30. 
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TABLE 3 

Aperture Correction Factors 



Galaxy Ha 3.6/im 4.5/^m 5.8^m S.O/im 24/im 



NGC 0024 


1, 


.03 


1, 


,02 


1, 


,02 


1, 


,02 


1, 


,03 


1.64 


NGC 0337 


4, 


,46 


4, 


,94 


4, 


,94 


5, 


18 


5. 


.69 


5.85 


NGC 0628 


1, 


,02 


1, 


,01 


1, 


,04 


1, 


05 


1, 


07 


1.64 


NGC 0925 


1, 


,09 


1, 


,07 


1, 


,07 


1, 


.07 


1, 


.08 


1.72 


NGC 1097 


2, 


,56 


2, 


,95 


2, 


,95 


3, 


,08 


3. 


.37 


3.21 


NGC 1291 


1, 


,16 


1, 


,23 


1, 


,23 


1, 


.25 


1, 


.31 


1.88 


NGC 1316 


4, 


,21 


4, 


,94 


4, 


,94 


5, 


.18 


5, 


.69 


5.85 


NGC 1512 


1, 


,47 


1, 


,61 


1, 


,61 


1, 


67 


1, 


78 


2.21 


NGC 1566 


2, 


,61 


2, 


,95 


2, 


,95 


3, 


.08 


3, 


.37 


3.21 


NGC 1705 


1, 


,00 


1, 


,00 


1, 


,00 


1, 


.00 


1, 


.00 


1.34 


NGC 2403 


1, 


,00 


1, 


,00 


1, 


,00 


1, 


.00 


1, 


.00 


1.34 


NGC 2841 


1, 


,55 


1, 


,61 


1, 


,61 


1, 


.67 


1, 


.78 


2.21 


NGC 2976 


1, 


,00 


1, 


,00 


1, 


,00 


1. 


.00 


1. 


.00 


1.34 


NGC 3031 


1 


,00 


1 


,00 


1 


,00 


1. 


.00 


1. 


.00 


1.34 


NGC 3184 


1, 


,16 


1, 


,23 


1, 


,23 


1, 


.25 


1, 


.31 


1.88 


NGC 3190 


4, 


,41 


4, 


,94 


4, 


,94 


5, 


.18 


5, 


.69 


5.85 


NGC 3198 


1, 


,47 


1, 


,61 


1, 


,61 


1. 


.67 


1. 


.78 


2.21 


IC 2574 


1 


,00 


1 


,00 


1 


,00 


1, 


.00 


1, 


.00 


1.34 


NGC 3351 


1, 


,10 


1, 


,15 


1, 


,15 


1, 


17 


1, 


21 


1.72 


NGC 3521 


1 


,16 


1 


,23 


1, 


,23 


1, 


25 


1, 


31 


1.88 


NGC '^627 


1, 


,15 


1, 


,23 


1, 


,23 


1, 


.25 


1, 


.31 


1.88 


NGC 3938 


1, 


,49 


1, 


,61 


1, 


,61 


1, 


.67 


1, 


.78 


2.21 


NGC 4254 


2, 


,60 


2, 


,95 


2, 


,95 


3, 


08 


3, 


37 


3.21 


NGC 4321 


1, 


,47 


1, 


,61 


1, 


,61 


1, 


67 


1, 


78 


2.21 


NGC 4450 


2, 


,20 


2, 


,11 


2, 


,10 


2, 


15 


2, 


22 


3.21 


NGC 4536 


1, 


,48 


1, 


,61 


1, 


,61 


1, 


.67 


1, 


.78 


2.21 


NGC 4631 


1, 


,08 


1, 


,07 


1, 


,07 


1, 


.07 


1, 


.08 


1.72 


NGC 4725 


1, 


,49 


1, 


,61 


1, 


,61 


1, 


.67 


1, 


.78 


2.21 


NGC 5194 


1, 


,02 


1, 


,04 


1, 


,04 


1, 


.05 


1, 


.07 


1.64 


Tololo 89 


2, 


,63 


2, 


,95 


2, 


,95 


3, 


.08 


3, 


.37 


3.21 


NGC 5474 


1, 


,01 


1, 


,01 


1, 


,01 


1, 


.01 


1, 


.01 


1.55 


NGC 5713 


4, 


,26 


4, 


,94 


4, 


,94 


5. 


.18 


5. 


.69 


5.85 


IC 4710 


1, 


,11 


1, 


,15 


1, 


,15 


1. 


.17 


1, 


.21 


1.72 


NGC 6822 


1, 


,00 


1, 


,00 


1, 


,00 


1. 


.00 


1. 


.00 


1.34 


NGC 6946 


1. 


,01 


1, 


,01 


1, 


,01 


1. 


.01 


1. 


.01 


1.55 


NGC 7331 


1. 


,49 


1, 


,61 


1, 


,61 


1, 


.67 


1. 


.78 


2.21 


NGC 7552 


4. 


,28 


4, 


,94 


4, 


,94 


5. 


.18 


5. 


.69 


5.85 


NGC 7793 


1. 


,00 


1, 


,00 


1, 


,00 


1. 


.00 


1. 


.00 


1.34 



Note. — Raw flux measurements are multiplied by these aper- 
ture correction factors to produce the corrected values. 
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Fig. 1. — 24 /im image of NGC 5194 (M51) with the final region selection indicated as small black circles, each 500 pc in diameter, along 
with several examples of background annuli (white circles) , each with an inner diameter of 500 pc and an outer diameter of 2000 pc. 
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38 39 40 41 38 39 40 41 

Log[i/L^(24)] or Log[L(Ha)] (CalzettI et oL 2005) 

Fig. 2. — Measurements of region luniinosities - fL,y(24) (eircles) and L(Ha) (diamonds) - for the four different tecfiniques deseribcd in 
the text compared with those of Calzetti et al. (2005) for a sample of matched 13 arcsecond (~517 pc) regions in NGC 5194. Final flux 
measurements were made using the mean of the global, local, and profile methods, as described in Section 4. 
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Fig. 3. — Histograms of ho^iv (24) / L{H a)) or Ajj^ for 1772 regions selected from the sample of 38 SINGS galaxies. The open 
histogram represents the entire region sample. The dot-dashed line is the fiducial cut proposed in the text at hog{vLi,(2A) / L{Ha)) = 2.7 
(A/fa ~ 3.3, Av ~ 4), corresponding to ~2 sigma above the mean of a Gaussian fit to the full distribution. The panel (b) is a cumulative 
histogram of the entire region sample; the dashed line is a fit to Log(i/L„(24)/L(//a))=[1.9,3.4] with high significance {B? = 0.98) yielding 
a power-law index of 1.4. There is no large sub-population of uniformly highly obscured regions in SINGS galaxies. 
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AfHa) 

2 3 4-557 




AfHa) 

2 3 4-557 



1 2 3 

Log[i/L,(24)/L(Ha)] 




1 2 3 

Log[t.L,(24)/L(Ha)] 



Fig. 4. — Attenuation distributions (Log{i'L,^(24)/ L{Ha)) or Aua) are shown separated by host galaxy 24 /im luminosity for the low 
inclination galaxies (i > 60°; (a)) and high inclination galaxies (i > 60°; (b)). In each panel the narrow-hatched histogram is the low 
luminosity bin {vLi,{24) < 8.4 X 10^^ [erg s^^]), while the wide-hatched histogram the high luminosity bin {vLu{2A) > 8.4 X lO*"^ [erg 
s~^]), and the open histogram includes all luminosities. The region attenuation distribution is more strongly peaJced for low inclination 
galaxies, while high inclination cases populate the wings of the distribution. The dashed line in each panel represents the median value of 
the distribution. Highly-obscured regions are found in both high and low luminosity galaxies. 
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Fig. 5. — Median attenuation for a representative set of galaxies as a function of aperture diameter used. The error bars indicate 
the standard deviation of all the region measurements included in the median. The dashed line indicates the size of the FWHM of the 
24 fj,ra PSF, below which the measurement is unresolved. The median attenuation is nearly constant with aperture size between 100 and 
1000 pc scales. In the sole case where we can probe below 100 pc, the dwarf galaxy NGC 6822, the measured attenuation appears to climb 
substantially; however, at these scales we are affected by the individual region morphology. For large apertures (>100 pc) the measured 
value of Log(i/Li/(24)/L(Ha)) is relatively insensitive to aperture size. 
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Log[F(24)/F(4,5)] 

Fig. 6.— (a) F(8.0)corr/F(4.5) versus F(24)/F(4.5) plot for the entire region sample, (b) F(8.0)corr/F(4.5) versus F(24)/F(4.5) plot for 
the highly-otscured subset, as defined by the attenuation cut described in the text. The contour encloses 95% of the total region sample. 
F(8.0) corr values have been corrected for the contribution of stellar emission, as discussed in the text. Flux densities are F^. quantities in 
erg Hz~^. Highly-obscured regions do not stand out from the rest of the region sample in terms of their infrared colors. 
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Fig. 7. — hog[uL^{24) / L(lia)] versus Log[;/L^ (24)] for the region sample, divided by host galaxy. The contour is common to all panels 
and encloses 95% of the total region sample, while the small points arc regions within the individual galaxy. The dotted line is the 
attenuation cut defined in the text, and the filled circles represent those regions in the galaxy that have been classified as highly obscured. 
The galaxy inclination is given in each panel. Highly attenuated regions tend to be more luminous. 
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Fig. 8. — Same as Figure 7. 
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Fig. 9. — Same as Figure 7. 
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Fig. 10. — Same as Figure 7. 
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Fig. 11. — Histogram of highly-obscured region locations in the host galaxy, given as a deprojected radial position in units of R2B 
(hatched histogram; multiplied by a factor of 10 for clarity). The open histogram shows the corresponding locations for the entire region 
sample. The dashed line shows the median location for the entire sample (0.38 R25) while the dotted line is the median location for more 
highly-obscured regions (0.26 R2b). More highly-obscured regions are generally located closer to the center of the host galaxy. 
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Fig. 13. — Same as Figure 12. 
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Fig. 14. — Same as Figure 12. 
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Fig. 15. — Dcprojcctcd radial profiles of Log(j^L^(24)/L(HQ)) for galaxies in the sample with inclinations greater than 60°. In cases 
where the region profile can be well fit by linear regression, the dashed line gives the linear fit to the radial range 0.0-0.65 R25- The dotted 
line corresponds to the attenuation cut (A^„ >3.3) applied in Section 6.1. 
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Fig. 16. — Same as Figure 15. 
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Fig. 17. — Zeropoint and slope measurements derived from radial trends in star-forming region attenuation versus the total galaxy 
H-band luminosity and galaxy inclination. The open squares correspond to low inclination (i < 60°) galaxies, while the filled circles are 
more highly-inclined cases. H-band luminosity errors are negligible; inclination errors are described in Section 2. Slope and zeropoint error 
bars are formal fitting errors. No strong correlations are seen; a tentative trend appears to exist between galaxy //-band luminosity and 
attenuation zeropoint and is corroborated by a similar trend between Mg and zeropoint. Inclination does not play a dominant role in 
governing the observed slope and zeropoint of the attenuation profile. 
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Fig. 18. — (a) Radial profile of Log(f L,y(24)/L(Ha)) for the entire region sample. While Zaritsky et al. (1994) Figure 7c showed little 
indication of a radial trend in Ay, with our larger sample a radial trend is evident though much less tight than in any individual galaxy, 
(b) Radial profile of Log(fLjy(24)/L(Ha)) for subset of our sample chosen to roughly replicate that of Zaritsky et al. (1994) with only a 
tentative indication of a radial trend. The apparent la<;k of a radial attenuation gradient in Zaritsky et al. (1994) was simply a consequence 
of plotting a region sample drawn from many galaxies, each with its own attenuation zeropoint. 



